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We recall the successes of the Hubbard-Stratonovich Transformation (HST) of many-body theory, 
point out its failure to cope with competing channels of collective phenomena and show how to 
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quantum field as suggested by the HST. 
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1. The Hubbard-Stratonovich transformation (HST) 
has a well-established place in many-body theory |l| and 
elementary particle physics Q. It has led to a good un- 
derstanding of important collective physical phenomena 
such as superconductivity, superfluidity of He 3 , plasma 
and other charge-density waves, pion physics and chiral 
symmetry breaking in quark theories [3(, etc. It has put 
heuristic calculations such as the Gorkov's derivation Q 
of the Ginzburg-Landau equations [H[ on a solid theoret- 
ical ground (fjj. In addition, it is in spirit close Q to the 
famous density functional theory [§] via the celebrated 
Hohenberg-Kohn and Kohn-Sham theorems Q. 

The transformation is cherished by theoreticians since 
it allows them to re-express a four-particle interaction 
exactly in terms of a collective field variable whose fluc- 
tations can in principle be described by higher loop di- 
agrams. The only bitter pill is that any approximate 
treatment of a many-body system can describe interest- 
ing physics only if calcuations may be restricted to a few 
low-order diagrams. This is precisely the point where the 
HST fails. 

Trouble arises in all those many-body systems in which 
different collective effects compete with similar strenghts. 
Historically, an important example is the fermionic su- 
pcrfhhd He 3 . While BCS superconductivity was de- 
scribed easily via the HST by transforming the four- 
electron interaction to a field theory of Cooper pairs, 
this approach did initially not succeed in a liquid of He 3 
atoms. Due to the strongly repulsive core of an atom, 
the forces in the attractive p-wave are not sufficient to 
bind the Cooper pairs. Only after taking the help of an- 
other collective field that arises in the competing para- 
magnon channel into account, could the formation of 
weakly bound Cooper pairs be explained [Ioj |. 

It is the purpose of this note to point out how to cir- 
cumvent the fatal fucussing of the HST upon a single 
channel and to show how this can be avoided in a way 
that takes several competing channels into account to 
each order in perturation theory. 

2. The problem of channel selection of the HST was 
emphasized in the context of quark theories in Q and 



in many-body systems such as He 3 in [fj. Let us briefly 
recall how it appears. Let x — (t, x) be the time and 
space coordinates, and consider the action A = Aq + A- m t 
of a nonrelativistic many-fcrmion system 



A-- 



r x [id t ~ e(-*v)] 



(1) 



where we have written ip x instead of i/j(x), and J for 
/ d A x, to save space. The symbol £(p) = e(p) — /tx de- 
notes the single-particle energies minus chemical poten- 
tial. Adding to A also a source term A s = J d 4 x(ip*r] x + 
c.c.) to form A = A + A s , the grand-canonical gener- 
ating functional of all fermionic Green functions reads 

The HST enters the arena by rewriting the interaction 
part with the help of an auxiliary complex field A X)X ' as 

Z[ri, r)*\ = J V<P*Vi(jVA*VA e ^[r ,4>A* A] +1 a s ( 2 ) 
with an auxiliary action 



vj* [id t - £(-»V)] 5 x ,x>^ 



-k&* x ,x'i>xi>x- - ilPx^x' &x,x' + HA^'lV^x' j, (3) 

Indeed, if the field A XjX ' is integrated out in ([5]), one re- 
covers the original generating functional. At the classical 
level, the field A XtX ' is nothing but a convenient abbre- 
viation for the composite pair field V x . x iij) x vb x i upon ex- 
tremizing the new action with respect to 5 A* x , , yielding 
5A/5A* X X , = (A XiX - - V XtX >ip x ip x >) /2V X , X > = 6. Quantum 
mechanically, there are Gaussian fluctuations around this 
solution which are discussed in detail in 0, Q . 

Expression ([3]) is quadratic in the fundamental fields 
ip x and reads in functional matrix form if x A xx >f X ' with 



[idt - CHV)] 5 X , 



A* , 

X.X 



[id t +aiV)]S Xy 



(4) 
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where f x denotes the fundamental field doublet ( "Nambu 
spinor") with f x = (ip x ,ip*), and f< = f* T , as usual. 
Since /* is not independent of f x , we can integrate out 
the Fermi fields and find 

Z[?f,r 1 }= J M*Me a[A *' A1 -54.'^[ G <»'V, (5) 

where j x collects the external source r\ x and its complex 
conjugate, j x = (r/ x , rj x ), and the collective action reads 



■A[A*,A] = -jTHog[iGX 1 ]+|/ |A X , X , 

J x.x' 



2 /V x , x >. (6) 



The 2x2 matrix Ga denotes the propagator iA 1 which 
satisfies the functional matrix equation 



[id t - e(-»v)] ^ 

- A* , 



^x,x' 



iidt + ti(iV)L x , 



[Ga]x',x" 



Gp Ga 
G t G, 



is equal to i6 x/x >>. Writing Ga as a matrix 

" r A _ U P 

the mean-field equations associated with this action are 
precisely the equations used by Gorkov Q to study the 
behavior of type II superconductors. 

With Z[r]*,r)] being the full partition function of the 
system, the fluctuations of the collective field A x>x / can 
now be incorporated, at least in principle, thereby yield- 
ing corrections to these equations. 

3. The basic weakness of the HST lies in the ambigu- 
ity of the decomposition of the quadratic decomposition 
of the interaction in ([1}. For instance, there exists an 
alternative elimination of the two-body interaction us- 
ing an auxiliary real field tp x , and writing the partition 
function as 



Z[rf,t)] = J Vi)*VipDipex-p[iA[i)*,i(),ip} + iA s ], (7) 
rather than ([2]), where the action is now 



A[i/J*,ip,<p}= \ip x [id t -Z(-iV)-ip(x)}5 X)X ,ip x , 

J x.x' L 

The new collective quantum field ip x is directly related 
to the particle density. At the classical level, this 
is obtained from the field equation 8 A/ dtp x = (p x — 
J dx'V X:X 'ifj*,tp x ' = 0. For example, iiV XyX > represents the 
Coulomb interaction (5 tjt //|x — x'| in an electron gas, the 
field ip x describes the plasmon fluctuations in the gas. 

The trouble with the approach is that when introduc- 
ing a collective quantum field A IjX ' or (p x , the effects of 
the other is automatically included if we sum over all fluc- 
tuations. At first sight, this may appear as an advantage. 
Unfortunately, this is an illusion. Even the lowest-order 
fluctuation effect is extremely hard to calculate, already 



for the simplest models of quantum field theory such as 
the Gross-Neveu model, since the propagator of the col- 
lective quantum field is a very complicated object. So it 
is practically impossible to recover the effects from the 
loop calculations with these propagators. Thus the use 
of a collective quantum field theory must be abandoned 
whenever collective effects of the different channels are 
important. 

The cure of this problem comes from the development 
some time ago, in the treatment of path integrals of var- 
ious quantum mechanical systems 
culation of critical exponents in (f. 
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and in the cal- 
-field theories fl2| . 
of a technique called Variational Perturbation Theory 
(VPT) [l3j]. This is democratic in all competing chan- 
nels of collective phenomena. The important point is 
that it is based on the introduction of classical collective 
fields which no longer fluctuate, and thus avoid double- 
counting of diagrams of competing channels by quantum 
fluctuations. 

4. To be specific let us assume the fundamental inter- 
action to be of the local form 



Aloe _ 9 



(9) 



where the subscripts "f, J, indicate spin directions, and 
we have absorbed the spacetimc arguments x in the spin 
subscripts, for brevity. 

We now introduce auxiliary classical collective fields 
and replace the exponential of the action in the gener- 
ating functional Z[rj,rf] = J T>ip*T>ip e l-A identically by 


with the new interaction 



-^int ~ 'Ant + r, / fx M x f x — 



■^a^^P^a (11) 



We now define a new free action by the quadratic form 
^gew = M _ h ljT Mxfx = h f x A^,f x ,, where /J 

denotes the fundamental field doublet /J = (tp a ,ip^). 
Then we rewrite „4g cw in the 2x2 matrix form analogous 
to Q as A?r =^-H x flM x f x = kf x A%£f x ,, with 
the functional matrix A^ x , being now equal to 



[id t -£(-i'V)]5 a /3+p a i3 

Kb \idt 



A Q ^ 

-£(*v)]<v 



PaB 



(12) 



The physical properties of the theory associated with the 
action A\°^ + A s can now be derived as follows: first we 
calculate the generating functional of the new quadratic 
action Aq bw via the functional integral ZQ CW [r),r]*} — 
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J T>ijj* Vifj e lA ° . From its derivatives we find the new 
free propagators Ga and G p . To higher orders, we 
expand the exponential e lAint in a power scries and 
evaluate all expectation values (i™/n!)([Ant w ]")o ew using 
Wick's theorem as a sum of products of the free parti- 
cle propagators Ga and G p . The sum of all diagrams 
up to a certain order g N defines an effective collective 
action A^ s as a function of the collective classical fields 

Obviously, if the expansion is carried to infinite order, 
the result must be independent of the auxiliary collective 
fields since they were introduced and removed in (fTTj) 
without changing the theory. However, any calculation 
can only be carried up to a finite order, and that will 
depend on these fields. We therefore expect the best ap- 
proximation to arise from the extremum of the effective 
action 



Ul,ll2|,ll7|. 



The lowest-order effective collective action is obtained 
from the trace of the logarithm of the matrix (fT2")l : 



Al - 



|Trlog 



(13) 



The 2x2 matrix G AyP denotes the propagator «[Af 

To first order in perturbation theory we must calculate 
the expectation value (Ant) 01 the interaction (fT2|) . This 
is done with the help of the Wick contractions in the 
three channels, Hartree, Fock, and Bogoliubov: 

For this purpose we now introduce the expectation values 

A «/3 = gtylipp), A f3a EE gtypipa) = [A* ap ]*, (15) 



and rewrite (Ant) as 



Pip = IPPcc]*, 



(16) 



(Ant)= (l/g) / (A| t A it - p n p it + p n p u ) 

J X 

- (l/2. 9 ) f (A Pa A} a + A% fi Ap a + 2p a p Pa p).(17) 

J X 

Due to the locality of A Q( g the diagonal matrix elements 
vanish and A a p = c Q( gA, where c Q( g is i times the 
Pauli matrix trig. In the absence of a magnetic field, 
the expectation values p a p may have certain symmetries: 

p n = p = p u , p n = p it = 0, (18) 

so that ([T7| simplifies to 



(A iat ) = (l/g) (|A| 2 +p 2 )-(AA* + A*A + 2pp) .(19) 

J X 

The total first-order collective classical action A l A p is 
given by the sum A A ^A Ap + (Ant)- 



Now we observe that the functional derivatives of the 
zcroth-ordcr action Aa p are the free-field propagators 
Ga, and G p 

[G A ] Q /3 , -, A° A = [G p ] a p . (20) 

OPap 



-A°, 



SA al3 '■ op a /3 

Then we can extremize Aa p w hh respect to A and p, 
and find that, to this order, the field expectation values 
([T6"l) are given by the free-field propagators (|20|) at equal 
arguments: 

A x = g[G A ]x,x, Px = g[G p ] x , x . (21) 

Thus we see that at the extremum, the action Aa p ' s 
the same as the extremal action 

A 1 [A,p] = A [A,p] 



9 Jx 



P 2 )- 



(22) 



Note how the theory differs, at this level, from the col- 
lective quantum field theory derived via the HST. If 
we assume that p vanishes identically, the extremum 
of the one-loop action At[A,p] gives the same result 
as of the mean- field collective quantum field action ([5]), 
which reads for the present (5-function attraction A\ [A] = 
A [A] — | J x |A| 2 . On the other hand, if we extremize the 
action A A p at A = 0, we find the extremum from the ex- 
pression A [A, p] = A [A, p] — ~ J x p 2 . The extremum of 
the first-order collective classical action (f22|) agrees with 
the good-old Hartree-Fock-Bogolioubov theory. 

The essential difference between this and the new 
approach arises in two ways: 

• First when it is carried to higher orders. In the col- 
lective quantum field theory based on the HST the 
higher-order diagrams must be calculated with the 
help of the propagators of the collective field such 
as {A x A x i) . These are extremely complicated func- 
tions. For this reason, any loop diagram formed 
with them is practically impossible to integrate. In 
contrast to that, the higher-order diagrams in the 
present theory need to be calulated using only or- 
dinary particle propagators Ga and G p of Eq. ([20)) 
and the interaction (fT2|) . Even that becomes, of 
course, tedious for higher orders in g. At least, 
there is a simple rule to find the contributions of the 
quadratic terms 5 J f x M x f x in ([TTj). given the di- 
agrams without these terms. One calculates the di- 
agrams from only the four-particle interaction, and 
collects the contributions up to order g N in an effec- 
tive action Av[A, p]. Then one replaces An[A,p] 
by Av[A — eg A, p — egp] and re-expands everything 
in powers of g up to the order g N , forming a new 
serie s T^, ^ n q' AJA, p]. Finally one sets e equal to 
l/<7 [la] and obtains the desired collective classical 
action An [A, p] as an expansion extending (122ft : 

N 

A N [A,p]=J2MA, P ]-(l/g) / (|A| 2 + A(23) 
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Note that this action must merely be extremized. 
There are no more quantum fluctuations in the 
classical collective fields A, p. Thus, at the ex- 
tremum, the action (|23p is directly the grand- 
canonical potential. 
• The second essential difference with respect to the 
HST approach is that it is now possible to study a 
rich variety of possible competing collective fields 
without the danger of double-counting Feynman di- 
agrams. One simply generalizes the matrix M x 
subtracted from ^4j nt and added to *4i n t in (fTT|) in 
different ways. For instance, we may subtract and 
add a vector field yj^a a ybS a containing the Pauli 
matrices a a and study paramagnon fluctuations, 
thus generalizing the assumption (|18j) and allow- 
ing for a spontaneous magnetization in the ground 
state. Or one may do the same thing with a term 
a a \7 l ipAi a +c.c. in addition to the previous term, 
and derive the Ginzburg-Landau theory of super- 
fluid He 3 as in @ . 

An important property of the proposed procedure is 
that it yields good results in the limit of infinitely strong 
coupling. It was precisely this property which led to the 
successful calculation of critical exponents of all 4> 4 the- 
ories in the textbook [l2[ since critical phenomena arise 
in the limit in which the unrenormalized coupling con- 
stant goes to infinity 18[. This is in contrast to another 
possibility, in principle, of carrying the variational ap- 
proach to higher order via the so-called higher effective 
actions (l9| . There one extremizes the Legendre trans- 
forms of the generating functionals of bilocal correlation 
functions, which sums up all two-particle irreducible di- 
agrams. That does not give physically meaningful re- 
sults [2(| in the strong-coupling limit, even for simple 
quantum- mechanical models. 

6. The mother of this approach, Variational Pertur- 
bation Theory 
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is a systematic extension of a vari- 
ational method developed some years ago by Feynman 
and the author It converts divergent perturbation 

expansions of quantum mechanical systems into exponen- 
tially fast converging expansions for all coupling strength 
[l7j |. What we have shown here is that this powerful 
theory can easily be transferred to many-body theory, if 
we identfy a variety of relevant collective classical fields, 
rather than a fluctuating collective quantum field sug- 
gested by the HST. This allows us to go systematically 
beyond the standard Hartree-Fock-Bogoliubov approxi- 
mation. 
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